Which group velocity of light in a dispersive medium? 
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Abstract 

The interaction between a light pulse, traveling in air, and a generic linear, non-absorbing and 
dispersive structure is analyzed. It is shown that energy conservation imposes a constraint between 
the group velocities of the transmitted and reflected light pulses. It follows that the two fields 
propagate with group velocities depending on the dispersive properties of the environment (air) 
and on the transmission properties of the optical structure, and are one faster and the other slower 
than the incident field. In other words, the group velocity of a light pulse in a dispersive medium 
is reminiscent of previous interactions. One example is discussed in detail. 

PACS numbers: 42.25.-p, 78.20.Ci, ll.30.-j 
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The concept of group velocity of an electromagnetic wave packet, born in the early 19 th 

n 

century nas been thoroughly analyzed in many propagation regimes: for dispersive, non- 
dispersive and absorbing properties of the propagation environment 2]. Such a quantity is 
naturally linked to the propagation velocity of energy. In fact, if the optical field propagates 
in dispersive and non-absorbing media, the two concepts of group and energy velocity tend 
to coincide In recent years this topic has induced a good deal of studies mainly within 
the context of the superluminal behavior of light, where the role of several phenomena, like 
material gain, Kerr nonlinearities, birifrangence, has been examined in detail [3]. Apart 
from all these interesting results, in our opinion there still exists some simple, but intrigu- 
ing, aspect in this scenario that deserves to be studied. Actually, in the present Letter we 
consider a wave packet incident over a generic dispersive, non-absorbing optical structure 
like a dielectric slab, a Beam Splitter or a finite Photonic Band Gap (PBG) etc.. On general 
grounds, due to the interaction of the wave packet with this optical structure, we have two 
emerging fields: the transmitted and the reflected pulse (as schematically depicted in Fig{T]) . 
We tackle the problem of their group velocities. As above mentioned, we will only assume 
that the propagation environment is dispersive, like air is, and that the involved field fre- 
quencies are far from the absorbtion lines. It will be shown that, due to the combination of 
the principle of energy conservation, the dispersive properties of the propagation environ- 
ment and the transmission properties of the generic optical structure, differences among the 
group velocities of incident, transmitted and reflected fields arise. For the sake of simplicity, 
and without loss of generality, we will refer only to fields having the wave vector k lying in 
the x, y plane, as shown in FigJTJ We can expand the electric field E^\x, y, t) of the incident 
wave packet, in plane waves: 

E^ l \x,y,t) = — [ E {i \cu)exp[i(k {i) ■ v - uot)]duo (1) 
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where r = (x,y). In Eq{T] = k(aj)s^\ where = (cos^, — sin#) and 9 is the angle of 
incidence, and 

k{u) = n{u)u/c. (2) 

In Eqf21 n(oj) is the refractive index of the environment which in nonabsorbing media is a 
real function of a real variable. As usual, the input-output behavior of a generic structure, 
with respect to the incident field, is describable by means of reflection and transmission 
coefficients, r{u) and t{oj) respectively, connected with the (intensity) reflectivity and trans- 



missivity through the well known relationships R(u) 
such hypotheses [gj], 9] 



r(uj)\ 2 and T(u) = \t(uj)\ 2 . Under 



T{uj) + R{uj) 



|f(a;)| 2 + |r(u;)| 2 = l. 



(3) 



The transmitted and reflected beams are 
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2vr J y ' 

exp[zk^ • r — ut)]du 
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E^(x,y,t) 



2tt J y ' y ' 
exp[i(k ( - r ' ) ■ r + Lut)]du 
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x 



(5) 



where k(' r ) = k(u)s^ r \ ky' = fc(w)s® and = (— cos 9, — sin 9) and = (cos 9, — sin 9). 
Group velocities are defined as 



where a = {i,r,t} and k$ a —k(uQ ) is the amplitude of the wave vector corresponding to 
the frequency ui^f 1 at which the spectrum E^ a \oS) has a maximum in the amplitude. This 
definition is due to the well known stationary phase method As we are interested in 
highlighting the effect of Eqj3] on group velocities, it is necessary to make some further 
hypotheses about the spectrum of the incident field as well as on the behavior of T(uS) and 
R(u). In particular, to exclude interactions producing very distorded wave packets, for which 
the definition of group velocity would be questionable, we will assume that the spectrum of 
the incident pulse is narrower than the transmission or reflection bands. This guarantees 
that the spectral range of the incident pulse is sufficiently limited to include at most one 
maximum (minimum) of T(u) (R(u)) or vice versa. However, also under such an hypothesis 
there is some pathological case that must be excluded. This happens when the frequency 
spectrum of the incident pulse has its maximum exactly at the point where T(u) = 1 (or 
R(u) = 1). In such a case the reflected (transmitted) pulse would be, in practice, either 
completely absent or very distorted. Generally speaking, if the modulus of the incident field 
\E®(u)\ has a maximum (for the intensity of the beams we are obviously interested only 
in moduli) for a given value of uj (determinated by the condition d\E^ l '{oj)\/dw = 0), which 
we denote as uj^ \ the reflected and transmitted spectra will show a shift in their maxima. 
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They will be respectively determined by the following two equations 

*y i« m +^i*<v)i=q (7) 

^M| r(a ,)| + ^|^ H | = (8) 
au; duo 

All the three fields would have a maximum in uo^ if, and only if, in that point there is a 
maximum (or a minimum) for T(u) or i?(co>) [11]. In such a case the fields would have the 
same group velocity. In all other situations, i.e. when dT{uj)lduj\ w = — dR(uj)/(ko\ « ^ 
(where and oJq denote the solutions of Eqs. [7] and [8] respectively), we can evaluate the 
moduli for the new group velocities as 

where a Taylor expansion has been made. In Eqq. [9] and [10] we have obviously k® ^ k$\ 
The reflected and transmitted pulses cannot have same group velocities. This property orig- 
inates from the observation that, by virtue of the constraint described by the Eqj3j spectral 
changes shown by the two emerging fields are not independent. In fact, it is easily seen that 
every time T(u) grows, in an interval of u>, then in the same interval R(u) must necessarily 
decrease and vice versa. This implies that they must always have different velocities, both 
different from the group velocity of the incident pulse. The difference between group velocities 
reported in Eqq. 0(721 depends on the dispersive properties of the propagation environment 
(the second derivative term) and on the properties of the optical structure considered (the k's 
difference term). Of course in any non- dispersive environment, i.e. when the modulus of 
wave vector is given by k = nu/c with n constant (for instance n — 1 for the vacuum), the 
second derivative terms in Eqs. I91TT01 are null. Therefore in such a case the above discussed 
phenomenon is absent. 

The competing effect of the group velocity dispersion (Av g ) must also be taken into account: 



It holds for the spread of each single wave packet due to the effect of the dispersive medium. 



This phenomenon tends to destroy the packet by spreading it while propagating [15|] • Hence 



only structures where this effect is under control will be considered. 

As an illustration let us calculate the difference for the group velocities of the two emerging 
pulses in air after an interaction with a dielectric slab. We can evaluate the second derivative 
term in Eqs. I9"][TU1 with an accuracy of 1 ppm [lsj by using the most accurate recent mea- 



surement of the refractive index of air 

a 



12] which improves Cauchy's formula for the same 
refractive index (j), n(u) = 1 + A(l + Bco 2 ) with A = 28.79 x 1(T 5 and B = 1.6 x l(T 33 s 2 . 
As to the optical structure, we consider a homogeneous, non-dispersive, rectangular dielec- 
tric slab with refraction index n = 1.5 and thickness d = 0.3/im along the x axis. The 
incident field impinging on this structure is supposed to be a broad-spectrum gaussian light 
pulse (white light), centrered around 0.555/xm: E^(u) = E^ exp{— (u — 0Jq ^) 2 /2a 2 }, with 
cJq = 3.39397 x 10 15 rad/s and er = 4 x I0 14 rad/s. It is well known that, in this case, the 
reflectivity 7Z± for the interface between air/glass depends on the polarization features of 
the field under study as well as on the (mean) angle of incidence. For the case of orthogonal 
polarization (electric field perpendicular to the x, y plane) 7Z± = sin 2 ^ — 6 t )/ sin 2 (6>j + 9 t ). 
Here Qi and 8 t represent the incident and transmitted angle, respectively. For 9{ = 65° it 
follows from Snell's law, 6t — 37° and TZ± = 0.228. The transmission properties of the entire 
optical structure are given by 

T ^ = IMF 1 - \ ^ 
1 + Fan (a) 

where F = 4TZ±/ (1—1Z±) 2 , a = uond cos(6 t ) / c [7]. In Figj2] (subplot a) we show the frequency 
spectrum of the incident pulse together with the transmission and reflection spectra of the 
dielectric slab (which is, after all, a Fabry-Perot multiple beam resonator), while in subplot 
b we show the frequency spectra for incident, transmitted and reflected fields [l^]. It is 
manifest that the transmitted and reflected pulses are shifted, with respect to the incident 
one, but in two opposite directions. Indeed we obtain: 

v f - v f = 2.47 x 10 2 m/s 
V <P ~ v 9 ^ = —1-65 x 10 2 m/s 

Recalling that the speed of light is known to an accuracy of ±lm/s, it is evident that such 
an effect is observable. 

Summarizing, we investigated the role played by energy conservation on the group velocity 
of a light pulse, propagating in dispersive media, after an interaction with a generic linear, 
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non-absorbing and dispersive structure. We have shown that in general, transmitted and 
reflected fields emerge with different group velocities. In a forthcoming paper the extension 
of the present formalism to the case of absorptive media will be considered. 
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FIG. 1: An optical beam incident over a linear dispersive and finite optical structure (such as a 
Beam Splitter or a PBG or a Lens etc.) is represented. The transmitted and reflected beams are 
schematically depicted as two other smaller wave packets, traveling in different directions. 
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FIG. 2: Transmitted (dotted line) and reflected (dashed line) field spectra originating from a 
rectangularly symmetric slab of thickness 0.3/um and refractive index n = 1,5. The incident field 
(continuous line) has an angle of incidence 6i = 65° . 
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